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Abstract 

We directly fit the QCD dimensional transmutation parameter, A^jg, to 
experimental data on e + e~ jet observables, making use of next-to-leading 
order (NLO) perturbative calculations. In this procedure there is no need to 
mention, let alone to arbitrarily vary, the unphysical renormalisation scale 
fi, and one avoids the spurious and meaningless "theoretical error" associ- 
ated with standard a s determinations. PETRA, SLD, and LEP data are 
considered in the analysis. An attempt is made to estimate the importance 
of uncalculated next-NLO and higher order perturbative corrections, and 
power corrections, by studying the scatter in the values of A-^g obtained for 
different observables. 
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1 Introduction 



The consistent extraction of the QCD coupling strength a s from experimen- 
tal data on a wide range of strong interaction processes has provided com- 
pelling evidence for Quantum Chromodynamics as the underlying theory of 
this sector of the Standard Model, for a recent comprehensive review see JT|. 
Notwithstanding attempts to exploit Lattice Gauge Theory calculations [2- 
4] , the bulk of these determinations extract a s by fitting the measured data 
to fixed-order perturbative predictions supplemented by computer- modelled 
hadronisation corrections. The asymptotic freedom of QCD leads one to 
hope that at sufficiently high energy scales (small enough values of the cou- 
pling) these determinations will be reasonably accurate. In this paper we 
shall be concerned with jet observables such as jet rates, thrust distribu- 
tions,..., studied in e + e~ collisions. For these observables next-to-leading 
order (NLO) perturbative corrections have been calculated [5-8], and non- 
perturbative hadronisation (power) corrections can be modelled using a vari- 
ety of computer Monte Carlo techniques || || . A difficulty with fixed-order 
renormalisation group (RG)-improved perturbation theory is that NLO pre- 
dictions depend on the dimensionful renormalisation scale \x which arises in 
removing ultraviolet divergences from the calculation. Usually one chooses 
\i proportional to Q the physical energy scale of the process, for instance the 
cm. energy in e + e~ collisions, so that fi = xQ, with x an unspecified dimen- 
sionless constant. This is then varied over some range around /i = Q, x = 1, 
the so-called "physical scale", say between fi = 10Q and /i = 0.1Q. The 
resulting fits to data of a s (xQ) are then customarily converted to a s {Mz) 
using the two-loop evolution of the coupling to yield a central value based on 
x = 1, and a "theoretical error bar" based on the variation between x = 10 
and x = 0.1. Unfortunately, the range of x to be considered is completely ad 
hoc, and there is no reason why the central values obtained should reflect the 
actual value of a s (Mz)- According to [l| this arbitrariness leads to "infinite 
discomfort in finite order" . The cure for this discomfort is to recognise that 
the dimensionful scale [i has no physical relevance whatsoever. The fact that 
fixed-order RG-improved perturbative predictions depend on it is a result 



of the standard way in which RG-improvement is carried out [10, 11 1, with 
the perturbation series being truncated at NLO, and a choice of scale which 
is Q-dependent. One should rather keep ji independent of Q. The renor- 
malisation scheme-dependent coupling a s (/i) then has no dependence on Q, 
the Q-dependence of the observable coming entirely from the perturbative 
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coefficients which contain "unphysical" /i-dependent logarithms and "phys- 
ical" Q-dependent ultraviolet (UV) logarithms. Fixed-order truncation of 
the perturbation series is no longer adequate, in particular such truncations 
with fj, independent of Q do not satisfy asymptotic freedom. The RG self- 
consistency of perturbation theory allows one to identify an infinite subset 
of RG-predictable logarithms at any given order of Feynman diagram cal- 
culation. If these are resummed to all-orders the "unphysical" /i-dependent 
logarithms contained in a s (fi) cancel against those contained in the pertur- 
bative coefficients, and one obtains a /i- independent result, which correctly 
builds the leading Q-dependence of the observable. This procedure has been 
termed "Complete Renormalisation Group Improvement" (CORGI) [|l(], 11], 
and it has been shown how to extend the argument to problems involving 
factorisation scales in addition to renormalisation scales |IT| . At NLO this 
approach yields exactly the same result as the Effective Charge approach 
of Grunberg |12|, |13| , which corresponds to choosing fi with standard RG- 
improvement in such a way that all the "physical" UV logarithms of Q are 
resummed. In general with the standard RG-improvement an infinite subset 
of these logarithms is omitted, with the result that one obtains /^-dependent 
results, and does not reproduce the correct physical leading Q-dependence 
of the observable. 

In fact it is possible to motivate this approach more straightforwardly by 
showing how QCD observables may be directly related to the dimensional 
transmutation parameter of the theory. In Section 2 we shall show how the 
dimensional transmutation parameter arises on the grounds of generalised 
dimensional analysis, modifying the analysis of [[14]]. The form of this rela- 
tion at the NLO level is then completely equivalent to the CORGI approach 
outlined above. The advantage of this derivation is that mention of the 
renormalisation scale fi and the renormalised coupling a s (n) can be essen- 
tially avoided, and the physical irrelevance of these quantities is manifest. 
In contrast the fundamental importance of the parameter A is stressed. In 
Section 3 we shall define the QCD e + e~ jet observables with which we shall 
be concerned. Section 4 contains direct plots of A^ obtained bin-by-bin 
in jet resolution criterion, thrust,. ..,etc, from hadronisation corrected data, 
using the direct relation between A^g- and the data. We attempt to esti- 
mate the uncertainty in A^g due to uncalculated NNLO and higher-order 
corrections, and possible power corrections, by looking at the scatter in the 
extracted values between different observables. In Section 5 we study the 
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dependence of the thrust distribution on the cm. energy Q using data span- 
ning the PETRA-LEP 1-LEP 2 energy range, and we perform joint fits for 
NNLO perturbative corrections and power corrections. Section 6 contains a 
discussion and our Conclusions. 



2 Direct Relation between A and QCD Ob- 
servable s 

In this section we shall derive from basic considerations of generalised dimen- 
sional analysis how the dimensional transmutation parameter A arises, and 
how it may be directly related to the QCD observable. The derivation owes 



much to the discussion of and closely follows that of [|15|]. We suppose 
that we have a dimensionless generic QCD observable TZ(Q), dependent on 
the single dimensionful (energy) scale Q. Quark masses will be taken to be 
zero throughout our discussion, the extension to the massive case has been 
considered in Since 7Z(Q) is dimensionless, dimensional analysis clearly 
demands that 

n(Q) = .(£), (i) 

where A is a dimensionful scale , which will turn out to be related to the 
dimensional transmutation parameter. There is an extra trivial possibility 
that 1Z(Q) = C, where C is a dimensionless constant. That is, there is 
no energy dependence. This trivial Q-dependence would be the case if the 
bare coupling of QCD was finite, since the QCD Lagrangian (with massless 
quarks) contains no massive parameters. Of course, in fact, the bare coupling 
is infinite, and an infinite renormalisation must be performed, leading to a 
functional relation as in Eq.(l). An obvious proposal is to invert Eq.(l) to 
obtain 

A = $^(Q)), (2) 
where $ _1 is the inverse function. This is indeed the basic motivation for 



Grunberg's method of Effective Charges [13|. We shall obtain the form of 



$ 1 by starting from the form of the derivative of 7Z{Q) with respect to Q, 
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imposed by dimensional analysis. We must clearly have, 



dft(Q) = B(K(Q)) 
dQ Q 



(3) 



where B(1Z(Q)) is a dimensionless function of TZ . This can be rearranged 
to obtain, 

This is a separable first-order differential equation. In order to solve it we 
will need to impose a boundary condition, and to know something of the 
behaviour of B(7Z). We shall assume that TZ has the perturbative expansion, 

K = a(l + Y,r n a n \ , (5) 

\ n>0 / 

where a=a s (n)/7T is the RG-improved coupling. The form of expansion in 
Eq. (5) can always be arranged by suitably scaling the observable and raising 
to an appropriate power. The required boundary condition will be given by 
asymptotic freedom, that is 1Z(oo) = 0. Integrating Eq.(4) one then obtains 

m ^r- = / t^t + K ■ ( 6 ) 
A n J B(x) 

The constant of integration has been split into lnA^ + k, where A-ji is a finite 
dimensionful scale specific to the observable TZ, and k is a universal infinite 
constant needed to implement lZ(oa) = 0. To determine k we need to know 
the behaviour of B(x) around x = 0. Returning to the perturbative series of 
Eq.(5) we note that the coupling a satisfies the beta-function equation, 

yf^- = /3(a) = -ba 2 \ 1 + ca + \ ^c n a n I , (7) 



<9/i 



+ J2 c ^ n ) 

n>l ) 



where b = (33 - 2N f )/6 , and c = (153 - 19iV»/126, are the first two 
coefficients of the beta-function for SU(3) QCD with Nf active (massless) 
flavours of quark. They are universal, whereas the subsequent coefficients 
q, (i > 1) are scheme- dependent. In fact, as shown by Stevenson |L7|], the 
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non-universal beta-function coefficients can be used to label the renormali- 
sation scheme (RS), together with the renormalisation scale. If we set fx = Q 
in Eq.(5), and differentiate with respect to InQ term-by-term using the beta- 
function equation (7), we can obtain B(lZ(a)) as a power series in a, finally 
if we invert the series 71(a) in Eq.(5) to obtain a(7Z) as a power series in TZ, 
we can obtain B(7Z) as a power series in 7Z . One finally finds for the series 
expansion of B(x) around x = 0, 



B(x) = -bx 2 ^1 + cx + ^y n x n j . (8) 

The first two coefficients b, c, are the universal beta-function coefficients. 
The higher terms Pi,(i > 1) are renormalisation scheme (RS)-invariant, and 



Q-independent, combinations of the and q. The first two are [13, 15, 17 



P2 = c 2 + r 2 - ric - r l 

Pz = c 3 + 2r 3 - 4r x r 2 - 2r x p 2 - r\c + 2r\ . (9) 

Knowledge of p n requires a complete N n LO perturbative calculation, that 
is a calculation of the for i = 1, 2, . . ., n and the q for i = 2, 3, . . ., n , 
in some renormalisation scheme, for instance MS. The fact that Eq.(8) has 
the same form as the beta-function equation (7) follows from the fact that 
there exists an RS in which 7Z = a, i.e. = 0,i > 0, and in this scheme 
the non-universal beta-function coefficients are q = pi, i > 1. The existence 
of this Effective Charge scheme |13[] is underwritten by the algebraic steps 
above from which (8) can be directly derived. Armed with knowledge of the 
form of B(x) around x = we see that the infinite constant of integration k 
will be of the form 

f c dx 

K = -/„ WY (10) 

where K(x) must be such that the singularity of 1/B(x) at x = in (6) is 
cancelled. This implies from (8) that 

K(x) = -bx 2 (l + cx + A(z)), (11) 

where A(x) is only constrained by the requirement that A(x)/x 2 is finite as 
x— >0. Different choices of the upper limit of integration, C, and the func- 
tion A(x), can be absorbed into the dimensionful constant A-^.. Convenient 
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choices are C = oo and A(x) = 0. With these choices (6) can be re- written 

as, 



dx f n ^ , 



The first integral on the r.h.s. of (12) gives 



b 1 

+ 



B(x) x 2 (l + cx) 



(12) 



cR 



1 + cJl 

Denoting the second integral by G(1Z) we have 

bln-Q- = F(K) + G{K) . 



(13) 



(14) 



The desired inverse function $ 1 of (2) can then be obtained by exponenti- 
ating (14), which gives 



F(K(Q))g(Tl{Q)) 



A 



n 



Q 



where T is the universal function 

^)=e"W = e^ m (l + l/cTZ) c/b , 

and 



g{JZ)=e- G ^' b = exp 



K(Q) ( i i 

dr I I 

, \B{x) bx 2 {l + cx) 



(15) 



(16) 



(17) 



If only a NLO perturbative calculation has been completed then our state 
of knowledge of B(x) is B(x) = —bx 2 {l + ex) since the NNLO and higher 
RS invariants p 2 ,p 3) ... of (8) will be unknown. From (17) we then have 
G(JV) = 1. We finally need to relate the observable-dependent constant of 
integration which arose on integrating (4), to the universal dimensional 
transmutation constant which depends only on the subtraction procedure 
used to remove the ultraviolet divergences, Aj^g for instance. Fortunately it 
turns out that they can be related exactly given only a one-loop (NLO) per- 
turbative calculation of the observable. To see this we begin by noting that 
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on rearranging (14) and taking the limit as Q^oo , we obtain an operational 
definition of A-^, 

A^ = lim Qexp(-F(K(Q))/b) . (18) 

Q— >oo 

We have used the fact that G(0) = together with asymptotic freedom. If 
we denote by a(Q) the MS coupling with [i = Q we see that it will satisfy 
the beta-function equation (7), of the same form as (4) for 71, with /%js(a) 
replacing B(1Z). This may be integrated following the same steps as above. 
The constant of integration A^ will be replaced by A^jg, and the coefficients 
Pi by the MS beta-function coefficients c^ s . Again choosing C = oo and 
A(x) = 0, we arrive at 

A m = lim Qexp(-F(a(Q))/b) . (19) 

From the perturbative expansion of 1Z in (5) we will have 

K(Q) = a{Q) + r{a{Q)f + . . . , (20) 

where we have defined for convenience r=r]^ s (/i = Q), as the notation sug- 
gests r is Q-independent. It is then straightforward to show that as Q^oo 

F(K)nF(a) - r + . . . , (21) 

where the ellipsis denotes terms which vanish as Q—>oo. Inserting this result 
into (18), and comparing with (19), one finally finds 

A^ = e r / b A m , (22) 

for the promised exact relation between the observable-dependent and uni- 
versal A's. The tilde over A is to draw attention to the fact that the above 
choice of infinite integration constant k does not accord with the standard 
choice [flSH ) which is based on an expansion of a in inverse powers of ln(/x 2 /A 2 ). 
This definition corresponds to translating k by the finite shift cln(6/2c), so 
that the standard Aj^g is related to Aj^g by 

f2c\~ c,b 

A m=[j) ^m- (23) 
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Finally assembling all this we arrive at the desired relation between the uni- 
versal dimensional transmutation parameter and the QCD observable 
K, 

Ams = QF(7l(Q)MnQ))e- r/b (2c/b) c/b . (24) 

Notice that all dependence on the subtraction convention chosen to remove 
ultraviolet divergences resides in the single factor e~ r ^ b , the remainder of the 
expression being independent of this choice. This is equivalent to the obser- 



vation of Celmaster and Gonsalves |19| that A's with different subtraction 



conventions can be exactly related given a one-loop (NLO) calculation. 
As noted above if only a NLO calculation has been performed then the state 
of our knowledge of the function B(1Z) in (4) is B(TZ) = —bTZ 2 (l + clZ), and 
then from (17) G{1Z) = 1. So at NLO the best we can do in extracting Ay^ 
from the data is 

A^ = QT(n(Q))e- r/b (2c/b) c/b . (25) 

If two-loop (NNLO) and higher-order perturbative calculations are available 
then Q(TZ) will differ from unity by calculable corrections. One can expand 
Qijt) as a power series in 71, 

g(K) = i-^n + o(n 2 ) + ..., (26) 

b 

where P2 is the NNLO RS-invariant defined in (9). Alternatively G(7Z) can be 
expanded in the exponent as a power series in TZ by expanding the integrand 
in (17), to give 

g{K) = exp[(-p 2 ft + 0{7l 2 ) + . . .)/&] . (27) 

One could also evaluate the integral in (17) numerically with B(x) truncated, 
so that at NNLO for instance B(x) = —bx 2 (l + cx + P2X 2 ). 

Focussing now on the NLO case where Q = 1 we note that (25) can be 
inverted to give 

K{Q) = ~c[l + W(z(Q))] 

z{Q) = (28) 



9 



where W is the Lambert VT-function [^TJ defined implicitly by the equa- 
tion W \z)exp(W \z)) = z. To be consistent with asymptotic freedom it is 
actually the W-\ branch of the function which is required [P^j . Eq.(28) is 
equivalent to the two- loop MS coupling with scale fi = e~ r ^ b Q, and in this 
scheme = 0. This scheme is sometimes referred to as the "Fastest Apparent 
Convergence" (FAC) scheme [17], and is equivalent to Grunberg's Effective 



Charge approach at NLO [12], [y|. Crucially, we have derived (28) without 
having to argue for a specific choice of scale. Starting from the form of Q- 
dependence of 1Z implied by dimensional analysis in (3), we simply solved 
this differential equation applying asymptotic freedom as a boundary con- 
dition. To define the required infinite constant of integration we needed to 
know the series expansion of B(x) around x = 0, Eq.(8), whose form is com- 
pletely scheme-independent, and we thus arrived at Eq.(14). The constant 
of integration could then be exactly related to the universal dimensional 
transmutation parameter A^g associated with use of MS subtraction to re- 
move ultraviolet divergences, given a NLO calculation of r, as in Eq.(22). 
In all of this the renormalised coupling a only ever appeared in interme- 
diate steps, playing, as neatly expressed in "the role of a conjuror's 
handkerchief- now you see it, now you don't!". This, of course, begs the 
question as to what is special about the Effective Charge (FAC) scheme, and 
why other choices of scale /i do not provide equally valid predictions for 1Z. 
The key is to identify the way in which the Q-dependence of 7Z{Q) arises. In 
the construction above it is built automatically by integration of (3), but how 
does it arise from the perturbation series in Eq.(5) ? The crucial observation 
is that the perturbative coefficients contain ultraviolet logarithms of Q. 
To see this we can rearrange Eq.(22) to obtain 

r = blnS- - bln-Q- , (29) 
Ams ^ 

or for a general choice of (MS) scale \i , 

ri (u) = bin J- - bln-Q- . (30) 
Ams ^ 

Thus ri is a difference of a scheme-dependent logarithm involving \x and a 
"physical" scheme-independent ultraviolet logarithm involving Q. In RG- 
improvement as customarily applied one chooses = xQ, and so the renor- 
malised coupling a(xQ) is Q-dependent. The perturbative coefficients r^fi = 
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Q), however, are Q- independent. The consistency of perturbation theory 
means that higher coefficients q are i th -order polynomials in r\ with Q- 
independent, but scheme- dependent coefficients |f[C|, [TTjl . Rearrangement 
of Eq.(9), for instance, gives 

r 2 = rl + nc + (p 2 - c 2 ) 

r 3 = r? + \cr\ + (3p 2 - 2c 2 )n + (| - |) . (31) 

From Eq.(30) we see that ri(/x = xQ) is Q-independent, and it follows from 
(31) that r 2 , r 3 , . . ., are therefore Q-independent too. Thus the Q-dependence 
comes entirely from the renormalised coupling, and is hence dependent on 
the unphysical renormalisation scheme parameter x. In contrast the idea of 
Complete RG-improvement (CORGI) JlO, 11] is to keep p strictly indepen- 



dent of Q. In which case the Q-dependence is built entirely by the ultraviolet 
logarithms of Q contained in the perturbative coefficients. Standard NLO 
fixed-order perturbation theory is then manifestly inapplicable, since one has 

Z 1 71 Q \ i i ^2 



Tl NLO = a(ji) + blnJ— - &ln-p- (a(/i)) z . (32) 
V Ams 

With /i constant, asymptotic freedom only arises if all the RG-predictable 
UV logarithms are resummed to all-orders. Given only a NLO calculation the 
RS-invariants Pi,(i > 1) are unknown, and so the resummation of the RG- 
predictable UV logarithms corresponds to setting the pi to zero in Eq.(31). 
If we further set c = and q = 0, (i > 1) , to simplify the analysis then the 
all-orders sum of RG-predictable terms reduces to a geometric series, 

K = a^) (l + Y^WM?) ■ (33) 

V n>0 / 

With these simplifications we will have a(fi) = l/61n(/i/Aj^g), and using 
Eq.(30) and summing the geometric series one obtains, 



H(Q) = a{n)/ 



l-(Hn ^L-Hn^. )„(/,) 

A MS A ^ 



l/b\n(Q/A n ) , (34) 



in which the unphysical /i-dependence has cancelled between a(/i) and the p- 
dependent logarithms contained in r\(p). In the realistic case with nonzero c 
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and Cj the simple logarithm of Q/An is replaced by the Lambert ^-function 
of Eq.(28). The key point is that the all-orders CORGI improvement can 
be carried out with any choice of \i to yield a /i-independent result. One 
has therefore directly traded unphysical /i-dependence for the physical Q- 
dependence. 

To further emphasise the connection of the suggested direct extraction 
of Aj^g with the standard approach we can consider the following result for 
Ajfg(ri,TZ), which we define to be the value of A^g obtained by fitting a NLO 
perturbative calculation in a scheme corresponding to the NLO coefficient r 1; 
to the data TZ . Notice that r\ completely labels the scheme at NLO. We can 
directly convert r\ into the MS scale \i since from (29) and (30) we have 

r x = r + bhij- . (35) 
It is then straightforward to derive the result [TJ] 

where / (ri , TV) is given by 

/(f „^W- ^ 1+ f^ )^, (37) 

In the CORGI approach r\ — and we have /(0, TV) = 0, so that the value 
of Aj^g obtained is Aj^g /Q (TV) , as expected comparing (24) and (25). Thus 
to the extent that Q(TV)~1 we obtain the actual value of Aj^g . As we have 
argued the estimate G(TZ)~l is the best we can do given only a NLO calcu- 
lation since we are in complete ignorance of the deviations of Q from unity, 
which will depend on the NNLO RS-invariant P2 of Eq.(9). Another way of 
saying this is that at asymptotic values of Q Eq.(25) will hold, and that the 
deviation of G(TZ) from unity provides an operational definition of how far 
from asymptotia we are, at Q = Mz, say. The scatter of the A^ values 
for different observables obtained from Eq.(25) thus provides unambiguous 
information about the size of sub-asymptotic effects (uncalculated NNLO 
and higher perturbative corrections and power corrections) . Variation of the 
renormalisation scale taking \x = xQ with the "physical scale" x — 1 giv- 
ing a central value merely serves to confuse matters. For instance taking 
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r = 10 and 7Z = 0.05, values typical of jet observables at Q = Mz, we find 
exp[/(10, 0.05)/6] = 2.44, and so using the "physical scale" the value of Aj^g 
extracted will be 2.44A^jg /Q(1Z). This will accurately determine A^jg if it 
fortuitously happens that £(7?.)~2.44 . We have, of course, no reason to 
suppose that Q differs from unity to such a drastic extent, or correspond- 
ingly that the effect of uncalculated NNLO and higher-order perturbative 
corrections, and possible power corrections should be so large. Varying the 
scale simply introduces an extra known factor into the determination of A^jg , 
which, other things being equal, i.e. if C?~l, will give values very different from 
the true one. 

In all this discussion we have considered strictly massless quarks. In reality 
the dimensional transmutation parameter has a dependence on the number 
of active quark flavours, Nf, so really we have A^f . The NLO correction r 
and the universal beta-function coefficients b, c, in Eqs. (24,25) also depend on 
Nf. Transformation between A^ for different values of Nf can be effected 

J MS J 

using the standard apparatus of the decoupling theorem augmented with a 
matching condition [03]. The matching condition has now been computed to 
the three-loop level |24" |. For all of our fits in Section 4, Nf = 5 will be the 



(5) 

active number of flavours, and we shall be extracting A^L 

' to MS 



3 Definition of the Jet Observables 

We restrict ourselves throughout to infrared safe observables and make all 
definitions in the centre-of-mass frame, with all sums running over N final 
state particles. We begin by reviewing the various possibilities for cluster- 
ing particles to form jets. Given a particular jet measure, y^ the following 
algorithm is common to all, 

1. Define a resolution parameter, y cut . 

2. For every pair of hadrons, hi and hj, evaluate the jet measure, y^. 

3. If the smallest occurrence of this quantity is less than the resolution 
parameter (i.e. min(?/jj) < y cu t) combine the corresponding hadron 
momenta, pi and pj into that of a pseudo-particle, pij according to a 
recombination prescription. 
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4. Repeat steps 2-4 until all hadrons and pseudo-hadrons have jet mea- 
sures greater than the resolution parameter. What remains are then 
denoted jets. 

By introducing a jet resolution parameter, y cut , we have made our defini- 
tion of a jet intrinsically infrared safe. Increasing its value permits a greater 
number of clusterings and thus few jet events are identified. Likewise, de- 
creasing its value finally results in all final state hadrons being assigned to 
separate jets. Within the theoretical framework, such small values probe 
deeply into the infrared region and thus require a thorough treatment of 
hadronisation. Jet measures are typically normalised by the total visible 
energy of the hadronic event, E v - ls , to give a dimensionless quantity. For a 
description of the multitude of different algorithms with a discussion of their 



merits see ||25|| . For this analysis we shall restrict ourselves to the JADE, 
Durham and Geneva jet finding measures applied by the various experimen- 
tal collaborations. 



The first jet measure to be proposed was by the JADE collaboration p6 
and simply uses, 

2EiEi(l -cos Oij) Mf 
IK ~ E. 



y j = *_^AL_ w (38) 



where Ei denotes the energy of a hadron, I, in the centre-of-mass frame and 
9i m is the opening angle of the pair under consideration. In the massless limit 
this measure corresponds to their invariant mass, Mfj. Having defined the 
jet measure we are still at liberty to define the procedure for recombining 
two hadrons into a pseudo-hadron. There are four immediately obvious pos- 
sibilities, denoted the E, EO, P and PO schemes. In all cases the subscript 
k denotes the pseudo-particle created by particles i and j. In the E scheme, 
we combine two particles according to their four-momenta, 

Pk=Pi + Pj- (39) 

Energy and momentum are explicitly conserved in this scheme. In the so- 
called EO Scheme the three-momenta of the pseudo-particle is rescaled to 
give it zero invariant mass, 

E k = Ei + Ej, (40) 

p k = nrx^T(P< + Pi)- (41) 
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As a result the total momentum sum of the event is not conserved. Conversely 
in the P Scheme we may conserve the total momentum of the event at the 
expense of the total energy conservation using 

Pfc = Pi + Pj, (42) 
E k = | Pk |. (43) 

Lastly we introduce a variation of the P scheme, the PO Scheme, by altering 
the jet measure such that after recombination, the total visible energy is 
changed such that, 



E vis = J2 E k- (44) 



Unfortunately the JADE jet measure turns out to introduce spurious clus- 
terings in certain circumstances whereby a resultant jet is formed in a direc- 
tion lacking any approximately collinear initial hadrons. This translates into 
theoretical problems when attempting to perform large infrared logarithm 
resummations where these correlations spoil the property of exponentiation 



in the two-jet limit fl27| . A subsequent attempt to suppress artificial recombi- 
nations within the jet clustering and hence improve its theoretical properties 
was suggested by Dokshitzer et al. [28|, termed the Durham or /^algorithm. 



It uses the minimum relative transverse momenta of two hadrons in the small 
angle limit, 

d 2mm(Ef,E])(l-cos6 lJ ) 

y l3 = m • ( 45 ) 

vis 

This form of clustering reduces the number of spurious recombinations and 
permits a straightforward theoretical implementation. As such it has now 
become the standard algorithm in use. We use the E scheme recombination. 
Lastly we consider a variant termed the Geneva algorithm proposed by 



Bethke et al. |29j that also attempts to reduce the spurious mis-clusterings 



of the Jade algorithm using the measure, 

g = ^ EjEjQ. - cosjy 
yij 9 {Ei + Ej) 2 ' { ' 

In contrast to the previous two proposals, the Geneva algorithm does not 
depend on the energy of the event, and has a preference to combine soft 
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particles with hard ones. This in turn reduces the correlations between soft 
gluons when performing infrared logarithm resummations. We also use the 
E scheme for recombination. 

With the jet finding algorithms in place we may now determine the n-jet 
rates (R n {ycut)) by the fraction of events with n resultant jets after clustering. 
We may then define the jet transition parameters, y n „ +1 that corresponds to 
the value of y cut where an event changes from (n + l)-jet-like to n-jet-like. 



We now turn to Event Shape Variables. Many of these variables are 
related and can be broadly categorised as follows. Note that they will in 
general contain both three and four-jet-like quantities. Thrust (T) is defined 
by maximising the net longitudinal momentum of final state particles along 



the direction of a thrust axis 30, pi 



T = max E ^ P V nT ' , (47) 

where p« denotes the final state particle momenta, and denotes the unit 
vector in the direction of the thrust axis, to be determined by maximising 
the above quotient. Defining for convenience r=l — T , we find that r varies 
between zero, for two back-to-back final state partons, up to a maximum of 
t = \ for spherical (isotropic) events. For planar events with three final-state 
partons, one finds a maximum value of r = | corresponding to a "Mercedes 
Benz" configuration. Two further variants, thrust-major (T ma j) and thrust- 
minor (T min ) can be defined. In T ma j the thrust axis is replaced in Eq . fl4~7|) 
by n ma j, which maximises the sum of momenta transverse to the thrust axis. 
In T min it is replaced by an axis n min which is the vector cross product of n T 
and n ma j. One can then define the oblateness O by 



0=T mai - T min . (48) 

Events can also be divided into two hemispheres (H a , H&) by a plane 
perpendicular to nr. We may then calculate the normalised invariant mass 
of each hemisphere (x = a, b) |33| . 



(49) 



Ml 

EL 
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This permits the possibility of four obvious combinations giving rise to 



Pt = Pa + Pb, (50) 

Pd = \p a -Pb\, (51) 

p H = max(p a ,p 6 ), (52) 

p L = min(p a ,p b ), (53) 

which correspond to the sum of jet masses, the difference of jet masses, 
the heavy jet mass , and the light jet mass, respectively. To lowest order in 
perturbative QCD, and assuming massless quarks, thrust and heavy jet mass 



are related by r = pn f34 | 



Other variants on thrust and jet masses are the jet broadening measures 
proposed in ||35|| . In each of the above hemispheres a and h one forms a jet 
broadening, B, by summing over the particles in that hemisphere, 

Bx = IP. x n r | (M) 

Once again we may compose a range of variables by the combinations, 

(55) 
(56) 
(57) 
(58) 

to make the sum of hemisphere broadening, the difference of hemisphere 
broadenings, the wide hemisphere broadening and the narrow hemisphere 
broadening respectively. For two-parton final states Bt = By/ = Br> = 0, 
and to lowest order in perturbation theory Bt = By/ = Bd = \0 = |T max . 

We can also define the so-called C and D-parameters from the eigenvalues 
of the infra-red safe linear momentum tensor [^, p7H , 

e mn = ^pMM , (59) 

where is the m-th component of the three- momentum p^, with i summed 
over all final state particles. As defined the tensor has unit trace. The 



Bx 


= B a + B b , 


B D 


= \B a — B b \, 


By/ 


= m&x(B a , B b 


B N 


= mm(B a , B b ) 



17 



C-parameter is then defined in terms of the eigenvalues of the tensor 6^, 
Ai, A 2 , A 3 by, 

C = 3(A 1 A 2 + A 2 A 3 + A 3 A 1 ) . (60) 

C = for back-to-back two parton final states and C = 1 for spherical 
(isotropic) events. For planar three-parton final states 0<C<| where one of 
the eigenvalues is zero. For values greater than | requires at least four final 
state particles. 

The D-parameter is defined by the combination 

D = 27X^3, (61) 

and only becomes non-zero for non-planar events (i.e. four or greater final 
states particles). 



Rather than describe an event by a single variable, we may consider 
inclusive two-particle correlations. The energy-energy correlation (EEC) 
[ 38] , |39 , fill is the normalised energy-weighted cross section defined in terms 
of the angle, Xiji between two particles i and j in an event, 

, Av 

1 r x+ 2 f f 

EEC{X) s m; E ./ ^ E W*u - *«> d *' < 62 > 

events 2 i,j vls 

where the argument, \i is t ne opening angle to be studied for the correlations, 
Ax, is the angular bin width and iV is the number of events. The angle x 
can be varied in the range 0° < x < 180° where the central region (x ~ 90°) 
is governed by hard gluon emission and the extremities (x ~ 0° and 180°), 
corresponding to collinear and back-to-back configurations, are expected to 
be sensitive to hadronisation. We may further define the asymmetric energy- 
energy correlation (AEEC) to be 

AEEC( X ) = EEC(180° - X ) - EEC( X ), (63) 

where now x is within the range 0°-90°. 

A recent addition to the set of e + e~ jet observables is the jet cone energy 
fraction (JCEF) |^TJ. Here the energy within a conical shell of opening angle 
X about the thrust axis is integrated, 

jcef{x) = j^- E f + J* E wr^ ~ < 64 > 

iV ^ X events J X~^ i ^vis 
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where 



Xi 



arccos 



Pi ■ n r 

IpI 



(65) 



is the opening angle between a particle, i, and the thrust axis vector, n T 
defined to point from the heavy jet mass hemisphere to the light jet mass 
hemisphere. The angle, \i * s within the range 0° < x < 180° and thus the 
hard gluon emissions will feature when \ > 90°. 

All the observables previously defined are either three or four-jet like. 
Additionally, we shall also investigate the two-jet rate. In the subsequent 
analysis we utilise experimental data at the Z°-peak from the SLD |42| . 
ALEPH [031, DELPHI H§ H§, L3 H§ and OPAL H7B collaborations. 



4 Direct Extraction of A 



MS 



In order to perform the direct extraction of A^jg- from the data for 1Z , 
using Eq.(25), we shall need to recast the perturbative expansions for the 
two, three and four-jet observables to be considered, so that they have the 
dimensionless form assumed in Eq.(5). For the three-jet-like observables, 71%, 
we will in general have a perturbative expansion at NLO of the form 

^(W-V^+W^Y- (66) 



2tt yr/ \ 2tt 

The two-jet case is equivalent to the three-jet case with the substitution 
IZ3 — > (1 — 7?.2)- Correspondingly four-jet-like observables, will have the 
expansion 

n m=Ai {^fj + B M {^-)\ (67) 

where A n denotes the //-independent tree-level coefficient of an n-jet-like 
quantity and B n the NLO coefficient. We may then by simple algebraic 
manipulation rewrite these in terms of the required dimensionless quantity 

as 

K 3 (Q) = 2 ^ 3(g) = <d + EM( <dY (68) 
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and 



W*(0)-\/-^r- - — + ^4-( — ) • (69) 



7T 



We are now in a position to calculate A^jg from Eq.(25) by substituting 
the experimental values of TZ(Q) and the fundamental quantity, r^/i = Q) 
which can be read from Eqs. (|68|) and (|69|) . We may then apply this to every 
experimental bin, enabling a direct extraction of across the kinematic 
range of the variable. In all cases we use the Monte Carlo programs EERAD 
48| and EERAD2 H9[ to calculate the NLO perturbative coefficients for three 



and four jet quantities respectively. 

Before we attempt to extract a value for A-^g there are a number of impor- 
tant issues worth considering. Firstly we must remember that even though 
we have defined a set of observables that attempt to reflect the underlying 
behaviour of the QCD partons, the effects of hadronisation will always be 
present. In some observables this will be more pronounced in certain regions 
of phase space resulting in the perturbative prediction failing to provide a 
reliable description. A number of Monte Carlo programs exist || [| that at- 
tempt to model this behaviour and have proved very successful. In addition 
non-perturbative power corrections have been studied phenomenologically 
and have displayed very positive results too. 

Secondly, from a purely perturbative QCD perspective, semi-inclusive 
quantities suffer from large kinematic logarithms at the exclusive boundaries 
of phase space. This manifests itself in the T\ parameter, growing in mag- 
nitude typically like — In 2 A for a variable A that goes to zero in the two-jet 
configuration. This drags the value for Aj^g to zero regardless of the true 
value, and indicates a breakdown in the NLO approximation since higher- 
order terms will be enhanced by powers of logarithms requiring an all-orders 
resummation. Furthermore, at the opposite end of the kinematic range we 
typically encounter a similar problem due to an end point in phase space. 
These occur when a variable goes from being n-jet-like to (n + l)-jet-like. 
Examples of this are the 1-thrust at | and the C-parameter at |. Above 
these values, the three-jet configurations do not contribute, resulting in the 
tree level term vanishing. Clearly this now upsets our definition of r\ since 
it diverges in a direction governed by the relative sign difference between LO 
and NLO in this limit. These characteristics can be seen in Fig. [3] for the 
1-thrust variable and Fig. |^ for the thrust minor variable. A more sophis- 
ticated way of handling the end point problem would be to define a new T\ 
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(1-T) Parameter 



Figure 1: The T\ parameter as a function of 1-thrust. It is calculated ac- 
cording to the ALEPH experimental bin sizes with errors attributed to the 
Monte Carlo integration of the LO and NLO coefficients only. 



value according to the ratio of the order o? s coefficient to order a 2 s in the re- 
gion corresponding to non-zero four-jet configuration contributions and then 
smoothly interpolate a value across the threshold. 

These two difficulties must be taken into consideration when attempt- 
ing to extract a value for A^jg . For a number of the three-jet quantities, 
hadronisation corrected data is analysed, providing a means of reducing that 
uncertainty In these cases the Jetset 7.4 hadronisation model || was im- 
plemented (unless otherwise stated) using bin-by-bin correction factors with 
errors estimated via statistical uncertainty. These factors were calculated as 



specified in [42[ and [45 
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Figure 2: The r\ parameter as a function of thrust minor. It is calculated 
according to the DELPHI experimental bin sizes with errors attributed to 
Monte Carlo integration of the LO and NLO coefficients only. 



Having identified the possible difficulties arising we must specify a set 
of criteria to perform a direct extraction of A^jg . This should concentrate 
on a plateau region of A^g in the central region of the kinematic range of 
the variable. We have adopted the following procedure for specifying the fit 
range, 

1. Decrease the fit range from its maximum value such that all r 1 values 
lie within a variation of 20% from the flattest region. 

2. Decrease the range further (if necessary) to the region where hadroni- 
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sation corrections are less than 40%. 



3. If more than three points are present, perform a single parameter \ 
fit to a flat line to calculate a value for with error. 



2 



4. Rescale the error according to ^Jx 2 J^df to obtain a A^jg prediction 
for the given variable and collaboration. 

If at any stage there are fewer than four consecutive bins surviving, the 
jet variable is considered unsuitable for the analysis. There is, of course, no 
guarantee that a variable will have a flat plateau over which to perform the 
fit. It may be such that the kinematic boundary effects dominate over the 
complete range. In these cases we are forced to disregard the variable. 

We have chosen to use an r\ criterion to avoid the problem of large kine- 
matic logarithms spoiling the fixed order perturbation theory. The param- 
eter, ri, clearly indicates the region where these logarithms are dominating 
the series and hence the breakdown of the NLO approximation. The value of 
T\ does not indicate where hadronisation effects may be considerable though. 
Therefore in order to give a proper treatment of the variables, we should 
use experimental data that has been corrected for hadronisation effects. We 
attempt to include a reasonably flat region across r 1 by allowing a 20% devi- 
ation from flatness (with errors taken into account). Since the r\ parameter 
varies smoothly across the kinematic range, this criterion permits a good 
measure of flatness. The value of 20% is chosen to tolerate minor deviations 
in r\ in the vicinity of the end points and any statistical fluctuations from 
evaluation of the NLO coefficients which are typically small. The resulting 
fit range should be relatively insensitive to small variations in the permitted 
percentage deviation. 

If hadronisation corrected data is available, we have adopted the proce- 



dure presented in ]|2] for excluding any bins that suffer from greater than a 
40% correction. 

Finally, we adopt a minimum x 2 test for fitting a flat line to the data 
points. The initial error (induced by Ax 2 = 1 from minimum) associated 
with the fit is then scaled by a/ x 2 l^df for Ndf degrees of freedom as promoted 
in the review of particle physics [PJ . This provides a value of Aj^g- for each 
observable measured by each collaboration. We must then consider how to 
combine the values. 

In considering the forthcoming fits, we must be careful not to underes- 
timate the errors. Dealing with different experiments' measurements of the 



23 



same observable will obviously have strong correlations. Typically the great- 
est difference between data sets will be due to statistical errors especially in 
the cases without any hadronisation corrections being applied. A procedure 
has been put forward by Schmelling termed the method of correlated 



averages, to combine correlated data when the exact correlation matrix is 
unknown. In this case, it is suggested that the degree of correlation is set 
by the \ 2 l^df value of the data set. In this way we are able to combine any 
number of correlated data without an unnatural reduction in the error. Sim- 
ilarly when combining errors with a \ 2 ' l^df greater than one we adopt the 
standard technique of rescaling the error by \Jx 2 /Ndf to improve the error 
estimation according to the quality of the fit. 

The experimental data for the three-jet observables (uncorrected for hadro- 
nisation effects) is taken from (|2| for SLD data, ||3| for ALEPH data, |44j 



for DELPHI data, [g§] for L3 data and g7| for OPAL data. Additionally, 
hadronisation corrected data is applied where available. In Figs. 3-15 we 
give the Aj^g plots for a set of thirteen observables for which hadronisation 
corrected data is available, and for which the fits to constant A^jg satisfy the 
criteria outlined above. The fit ranges are indicated by vertical dashed lines, 
and the best overall fit to Aj^g for the data sets considered, by a horizontal 
dash-dotted line. These values of A^g are assembled in Fig. 17. 
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Figure 3: The 3 Jet Rate using OPAL data 
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Figure 4: The Wide Jet Broadening using DELPHI and SLD data 



400 



300 



200 



100 



The Total Jet Broadening 



+ Delphi (JETSET 7.4) 
X Sid (JETSET 7.3) 



■+-X 



_+_._.* JX-+-. 



0.05 



0.1 



0.15 



0.2 



0.25 0.3 
Bsum 



Figure 5: The Total Jet Broadening using DELPHI and SLD data 
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Figure 6: The Sum of Hemisphere Masses using DELPHI data 
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Figure 7: The Heavy Hemisphere Mass using DELPHI and SLD data 
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Figure 8: The C Parameter using DELPHI, OPAL and SLD data 
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Figure 9: The 1-Thrust Parameter using DELPHI, OPAL and SLD data 
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Figure 10: The 2^3 Jet Transition Parameter using DELPHI and SLD 

data 
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Figure 11: The 2-^3 Jet Transition Parameter using DELPHI data 
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Figure 12: The 2^3 Jet Transition Parameter using SLD data 
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Figure 13: The 2 — > 3 Jet Transition Parameter using SLD data 
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Figure 14: The 2^3 Jet Transition Parameter using SLD data 
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Figure 15: The Jet Cone Energy Fraction using DELPHI and SLD data 
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We note that some of the variables described in Sec. 3, in particular the 
difference in hemisphere masses, and the EEC and AEEC, have not been 
included, since the fitting criteria outlined earlier are not satisfied. For com- 
parison the corresponding fitted values using data not corrected for hadro- 
nisation effects is given in Fig. 16. Without hadronisation corrections the 
fitted central values of lie in the range between 50 MeV and 275 MeV. 
Including hadronisation corrections, reduces this range to between 100 MeV 
and 200 MeV. This factor of two uncertainty in A^jg is then to be taken 
as providing an estimate of the likely size of remaining uncalculated NNLO 
and higher-order perturbative corrections. If we assume that for some of 
the observables G(7Z)~1, then the factor of two uncertainty in A^jg can be 
translated into an estimate of the potential size of the two-loop NNLO RS- 
invariant p 2 . Using Eq.(17) for Q(JV) expanded as in Eq.(27), we see that 
G(7l)~2, with the typical value of 7?.~0.05, corresponds to p 2 ~50. Without 
explicit NNLO calculations there is no rational basis for assigning a central 
value for Aj^g with an error. What one can say, however, is that inclusion 
of hadronisation corrections does reduce the spread in the extracted A^g , 
and that this remaining spread indicates a value of p 2 which is not so large. 
We have no basis for estimating how big we might expect p 2 to be for these 
observables. We can note that for the e + e~ i?-ratio which has been computed 



to NNLO, p 2 = —15.1 for N* = 5 active flavours [PJ. We can convert A 



L MS 



to a s [Mz) using the two-loop beta-function equation. The corresponding 
values for hadronisation corrected data is given in Fig. 18. 
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Summary of Measurements 

Uncorrected for Hadronisation Effects 
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Figure 16: Summary of A^jg measurements uncorrected for hadronisation 
effects. 



32 



Summary of Measurements 



Corrected for Hadronisation Effects 
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Figure 17: Summary of Aj^g measurements corrected for hadronisation ef- 
fects. 
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Summary of a § Measurements 

Corrected for Hadronisation Effects 
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Figure 18: Summary of a s measurements corrected for hadronisation effects. 
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We finally present the Aj^g values obtained by fitting to the data on four- 
jet observables. We use the DELPHI data of in all cases. In this case we 
only analysed data uncorrected for hadronisation effects. The results for five 
variables are given in Fig. 19. With the exception of the Light Hemisphere 
Mass and the Narrow Jet Broadening the values are grouped between 150 
and 200 MeV, reasonably consistent with those obtained from the three-jet 
observables. It may be that the Light Hemisphere Mass and Narrow Jet 
Broadening have rather large hadronisation corrections. 

Summary of Measurements 

Uncorrected for Hadronisation Effects 
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Figure 19: Summary of A^jg measurements for four jet observables uncor- 
rected for hadronisation effects. 
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5 Energy Dependence of the Thrust Distribu- 
tion 



In this section we consider how to include power corrections in our formalism. 
It is widely accepted that physical observables in general will be subject to 
"non-perturbative" power-like corrections in the hard interaction scale, Q. 
That is to say, there will be terms contributing to cross-sections that cannot 
be expanded out in the typical perturbative manner arising from expressions 
of the form 

A\ k f-k 



QJ exp feJ' (70) 

Perturbative techniques cannot describe these terms accurately but have 
made attempts at predicting the leading behaviour to the power corrections 
via renormalon-inspired analysis |52|, [53[] and dispersive techniques []54] . Tak- 



ing the generic form of these power corrections, we can alter our perturbative 
expansion for a dimensionless observable 71, 

TZ = a + na 2 + r 2 a 3 + ... + ^(1 + A x a + X 2 a 2 + A 3 a 3 + ...), (71) 

where we have assumed a leading power correction with exponent 1 (i.e. ^). 
To include these term in the analysis of Sec. 2, we must take the deriva- 
tive with respect to Q. We may then rewrite Eq.(8) incorporating power 
corrections as 

AID 

B{n{Q)) = — — = -bn 2 (i + cn + p 2 n 2 + P3 n 3 + ...) 

dmQ 

+ KoTZ-^e-^il + k x K +...) + ... , (72) 

where the be related to the X n . For example, the leading power 

correction coefficient gives a — \q/Q contribution to the /3-function. Using 
Eq.(12) to get Q in the leading approximation of 1Z we find, 

Q w A n e 1/bn (c7Z) c/b . (73) 

Substituting this back in we obtain 
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where we have converted A R to Aj^g . Having made the connection between 
Ko and Ao we may incorporate the power correction term into Eq. (24) via the 
Q(7l(Q)) function given in Eq.(17). Expanding out Q(7l(Q)) to the accuracy 
of NNLO and leading power corrections gives 

g(K(Q)) ~ 1 - + -f-j^^- c/b - (75) 

Substituting this back into Eq.(24) we finally obtain 

Ams = QFWfi) (l - + T-W- 11 ^) ^ (f)^ • (76) 

We are fortunate enough to have experimental measurements of the thrust 
distribution at a wide variety of energy scales from PETRA to LEP2. We 
shall be considering the following data 

• PETRA (Detector- Tasso, Facility-DESY) |5§ 

— Centre of Mass Energies - 14, 22, 35 GeV 

• PEP (Detector- Mark-II, Facility-SLAC) @] 

— Centre of Mass Energy - 29 GeV 

• TRISTAN (Detector- Amy, Facility-KEK) g7| 

— Centre of Mass Energy - 52 GeV 

• SLC (Detector- SLD, Facility-SLAC) 

— Centre of Mass Energy - 91 GeV 

• LEP (Detectors- Aleph, Delphi, L3, Opal, Facility-CERN) 

— Centre of Mass Energy - 91 GeV 

Unfortunately as the LEP 2 data suffers from large errors due to poor statis- 
tics, we are forced to exclude it from the analysis. We would expect that 
the Ajyfg directly extracted from Eq.(25) should approach the actual value as 
the energy increases, and sub-asymptotic effects become smaller with Q(7Z) 
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Figure 20: The Energy Dependence of Aj^g Measurements from the Thrust 
Distribution 



approaching more closely its asymptotic value of unity. This trend can in- 
deed be seen from Figs. |20] and [21] (where we have removed the error bars for 
clarity) . In the limit 1 — T — > | we are faced once again with the problem of 
the kinematic end point dragging the value to zero. In the 1 — T — > limit, 
large kinematic logarithms dominate from the emission of soft and collinear 
gluons. Clearly in between we see a "flattening" of the value with 1 — T. 
This can be interpreted as higher order and power-like corrections having less 
influence, and hence the NLO approximation becoming more reliable. 

It is possible to consider the formalism from a different perspective whereby 
we accept that Aj^g is a constant and that any deviations from its true value 
are due to higher order effects and power corrections. In this case we may 
then try to approximate these contributions by the leading terms as given by 
Eq. (|76|). In this section we highlight a simple mechanism for investigating 
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these effects. We stress that this by no means provides a highly accurate 
estimate of these terms, but merely an indication of how important they are. 
We begin by rewriting Eq. (ff6|) as 



K = M(n,r)p 2 + C(n,r), (77) 

where M(lZ,r) and C(lZ,r) are known terms at NLO, having specified a 
value for A^jg . This is simply the equation of a straight line in (n , p 2 ) 
space. Since these two quantities are Q- independent, we may plot the lines 
corresponding to different centre-of-mass energies and expect them to cross 
over at the solution. This is illustrated in Fig. [22] where we have taken a value 



of AjyFg = 200MeV. We have made no attempt to incorporate the errors, 
which for the LEP 2 data will be considerable. In all cases, the central value 
is taken. This naive procedure does give a promising result, though. There 
appear to be two predominant localised cross-over regions where the lines 
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200MeV 



appear to converge. These correspond to p<i and n both small or p 2 ~~ 40 
and tvQr 



0.5. Similar fits for the < 1 — T > were performed in [[58 



6 Discussion and Conclusions 



In this paper we have shown that one can directly relate QCD observables to 
the underlying dimensional transmutation parameter of the theory, leading 
to the relation of Eq.(24). This derivation does not require the specification 
of a renormalisation scale, or use of the convention-dependent renormalised 
coupling a s (fi). It is obtained simply by integrating up the obvious dimen- 
sional analysis statement of Eq.(3). This integration corresponds to identi- 
fying and resumming the complete set of ultra-violet logarithms of Q which 
reside in the perturbative coefficients. In standard RG- improvement an in- 
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finite subset of these is tacitly omitted, and this results in the problematic 
scale-dependence of fixed-order predictions which plagues standard analyses. 
The only approximation in the present approach is the incomplete knowl- 
edge of the function Q(TZ) which approaches unity for asymptotic Q, and 
deviates from it by two-loop (NNLO) sub-asymptotic effects characterised 
by the NNLO RS-invariant p 2 , and by power corrections. Given that these 
are presently unknown for the QCD e + e~ jet observables we wished to anal- 
yse, the best we could do was to use Eq.(25), corresponding to the asymptotic 
expectation G(Jl(Q))~l. As we discussed in Sec. 4, to avoid regions contain- 
ing large kinematical logarithms arising in the two-jet limit, and kinematical 
endpoints, it was necessary to choose the range in the observable over which 
the fits to constant Aj^g were performed, rather carefully. Notwithstand- 
ing the fact that these fitting criteria are inevitably somewhat ad hoc, the 
procedure crucially avoids the additional scatter in extracted A^jg values re- 
sulting from the arbitrary variation of the renormalization scale, which is 
customarily employed. One finds that the range of directly fitted Aj^g val- 
ues obtained is significantly reduced using hadronisation corrected data (see 
Figs. 16 and 17). The remaining scatter in A^jg values displayed in Fig. 17 is 
then to be attributed to the presence of uncalculated NNLO and higher-order 
sub-asymptotic effects in Q(TZ). The factor of two scatter in Aj^g, roughly 
between 100 and 200 MeV, then corresponds on using the exponentiated 
form for G(TZ) in Eq.(27), to an estimate |p2|~50, indicating significant but 
not huge uncalculated two-loop effects. A firmer statement about the value 
of A-^g is only possible once these effects have been computed. There seems 
to us to be no point in extracting a purported central value of Aj^g together 
with an error in the meantime. 

In Sec. 5 we performed a complementary exercise in which, assuming a 
fixed value of A^jg = 200MeV, the Q-dependence of the thrust distribution 
over the energy range from PETRA to LEP 2 was used to fit for the NNLO 
RS-invariant p 2 and power corrections. As seen from Fig. 21 as Q increases 
the directly extracted A^g is reduced, and its distribution in thrust flat- 
tened, as would be expected if the sub-asymptotic effects contained in G(TZ) 
are becoming smaller. The fits for p 2 and the k,q, the parameter controlling 
1/Q power corrections, shown in Fig. 22 reveal two possibilities. One corre- 
sponds to both p 2 and kq small, and the other to k q small and p 2 ~ — 50, a 
value consistent with that estimated from the scatter of A^g values in Fig. 17. 
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There are two major directions in which the formalism described here 
could be improved. The first would be to resum to all-orders large kinematical 
infra-red logarithms in the function Q{TZ). This would allow fitting over a 
much larger range in the observables. Whilst in principle straightforward a 
number of technical complications present themselves. The second would be 
to generalize the direct extraction of A^g to see how Eq.(24) is modified if 
one has massive quarks. We hope to report developments in both these areas 
in future papers. 
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